Optimal control of a qubit coupled to a non-Markovian environment 
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A central challenge for implementing quantum computing in the solid state is decoupling the qubits 
from the intrinsic noise of the material. We investigate the implementation of quantum gates for a 
paradigmatic, non-Markovian model: A single qubit coupled to a two-level system that is exposed 
to a heat bath. We systematically search for optimal pulses using a generalization of the novel open 
systems Gradient Ascent Pulse Engineering (GRAPE) algorithm. We show and explain that next 
to the known optimal bias point of this model, there are optimal shapes which refocus unwanted 
terms in the Hamiltonian. We study the limitations of controls set by the decoherence properties. 
This can lead to a significant improvement of quantum operations in hostile environments. 



A promising class of candidates for the practical re- 
alization of scalable quantum computers are solid state 
quantum devices based on superconductors [U El [3j HI [5] 
and lateral quantum dots [B]. A key challenge to over- 
come in this enterprise is the decoherence induced by the 
coupling to the macroscopic bath of degrees of freedom 
not used for quantum computation (see e.g. ref. [7J for a 
recent review). Many of these decoherence sources can 
be engineered at the origin. In the case of intrinsic slow 
noise originating from two-level fiuctuators (TLFs) this is 
much harder [HUH], albeit not impossible [TO] EI]- Thus, 
in order to advance the limitations of coherent quan- 
tum manipulations in the solid state, it is imperative 
to find strategies which accomodate this kind of noise. 
A number of methods have been proposed by intuition 
and analogies to different areas, such as dynamical de- 
coupling O |T3], the optimum working point strategy 
[U [14] , and NMR-like approaches [15] . Even in fight 
of their success, it is by no means clear, whether even bet- 
ter strategies can be formulated and, on a more general 
level, where the limits of quantum control under hostile 
conditions are reached. 

We resort to numerical methods of optimal control. 
The closed systems GRAPE (gradient ascent pulse en- 
gineering) algorithm |16j has been proven useful in spin 
and pseudo-spin systems |17j . an important example of 
the latter being coupled Josephson devices [TS]. It was 
recently extended to open yet strictly Markovian systems 
[19] . Other recent optimal control results also include the 
presence of noise and decoherence. General pure dephas- 
ing was considered in [2D]. Ref. ET] treated a semiclas- 
sical random-telegraph noise (RTN) model in the high- 
temperature limit, while ref. [22 focussed on two qubits 
and classical 1/f noise. Optimal state transfer in the 
spin-boson model was considered in ref. [23] . Ref. [24] op- 
timized qubit gates in the presence of finite-dimensional, 
dissipation-free environments. In this Letter, we general- 



ize the GRAPE algorithm to include a complex environ- 
ment leading to non-Markovian qubit dynamics and non- 
Gaussian noise. We show that next to an optimal work- 
ing point there is also an optimal pulse shape and optimal 
gate duration. Accelerating the fluctuations can improve 
the gate fidelity. We discuss the physics ultimately limit- 
ing the gate performance no longer correctable by pulse 
shaping. 

Model and method.— In macroscopic samples, 1// noise 
[H [5J E5j occurs in most observables. The Dutta-Horn 
model EEJ EI] explains this phenomenon by the (clas- 
sical) superposition of TLFs which randomly jump be- 
tween their states, a process known as random tele- 
graph noise (RTN). In small, clean samples the discrete 
nature of the noise process from a single dominating 
fluctuator [TJ EE] can be resolved. This leads to semi- 
phenomenological Hamiltonians E3 [501 HI] ■ 

We specifically model a qubit coupled to a single TLF 
by H = Hs + Hj + Hb ■ Hs consists of the qubit and the 
coupled two-state system, i.e. 

H S = E 1 (t)<T z + A<7 x +E 2 T z +AcT z T z . (1) 

<7j and T, are the usual Pauli matrices operating in qubit 
and fluctuator Hilbert space respectively. Ei(t) is time- 
dependent and serves as an external control. The source 
of decoherence is the coupling of the fluctuator to the 
heat bath, which leads to incoherent transitions between 
the fluctuator eigenstates, 

H I = Y,^(T + h + T-b]), H B =^ThUiblbi. (2) 

i i 

We introduce an Ohmic bath spectrum J(oS) — 
J2i — w i) — Kuj<d(u> — uj c ) containing the couplings 

A,;, the dimensionless damping k, and a high-frequency 
cutoff lu c (which we assume to be the largest frequency 
in the system). 
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We are mostly interested in the qubit evolution in the 
limit of slow TLF flipping. As a first step we describe the 
dynamics of the larger qubit <g) TLF system by a master 
equation, tracing out the bath along the lines of [71 152], 
We arrive, in the motional narrowing regime fcgT > KE2, 
at the Bloch-Redfield equation 

P(t) = ^[Hs,p(t)} + [T + ,^p(t)} + [r-,j:+p(t)} 

-[t- p(t)E+]-[r+ p(t)£o] (3) 
with the different rate tensors (s = 0, 1) 



{ihf 



= '— I -//' / dujJ(u){n(uj) + s)e ±iut ^(t'). 

(4) 

Here, n(o>) is the Bose function. Note, that the rate 
tensors explicitly depend on the control E\ (t) due to the 
interaction representation of the operators in eq. [4] 
Since tracing out the TLF at this stage would lead to an 
intricate non-Markovian master equation, we treat the 
qubit-TLF interactions exactly, i.e. the rate tensors act 
on the combined qubit-TLF system. 

Our model goes far beyond a simple RTN noise model 
[13j and captures the correlations between qubit and TLF 
[29l [31] . Still, it is useful to introduce the parameters 
of the RTN which would result for A -> 0. The TLF 
flipping rate is 7 = 2kE 2 coth(E 2 /T), the sum of the 
excitation and relaxation rate. It enters the two-point 
noise spectrum of random telegraph noise 



dt e 



(r z (t)T z (0)) 



cq 



A 2 



7 



uj + 7 



(5) 



This is the Fourier transform of the interaction represen- 
tation of t z assuming the bath in equilibrium. In this 
limit, we can find relaxation rates 1/T± = ^S{2E) and 

1/T 2 = l/2Ti + f|s(0) with£= ^A 2 + Ej. It has been 
shown in jS] that pure dephasing can be suppressed by 
keeping \E\ \ <C A during all manipulations, the optimum 
working point strategy. Additionally in the slow flipping 
regime, l/Tj-relaxation is small since S(2A) ~ A 2 7/4A 2 . 

We formulate the control approach by rewriting 
the master equation ^ as p(t) = — (iH(Ex(t)) + 
T(E\(t)))p(t) with the Hamiltonian commutator super- 
operator H(Ei(t))(-) = [H(Ei(t)),-] and the relaxation 
superoperator T, both time-dependent via the control 
Ei(t). The formal solution to the master equation is a 
linear quantum map operating on a physical initial state 
according to pit) — F(t)p(0). Thus F itself follows the 
operator equation of motion 



f = - [m + r) f 



(6) 



with initial condition F(0) = 1, as in rcf. 19]. Here, 
multiplication of quantum maps denotes their concate- 
nation. The task is to find control amplitudes E\(t) with 



t £ [0, t g ], t g being a fixed final time, such that the dif- 
ference SF = Fjj — F(t g ) between dissipative time evo- 
lution F(t g ) obeying eqn. [6] and a target unitary map 
Fjj is minimized with respect to the Euclidean distance 
||<5F|| 2 = tr {SF^SF}. Clearly, this is the case, when the 
trace fidelity 



Retr 



(4F(t g )} 



(7) 



is maximal. Note, that in an open system, one cannot 
expect to achieve zero distance to a unitary evolution 
Fjj [15] . The goal is to come as close as possible. 

We find the pulses by a gradient search. We digitize 



nt g ) 



N 



■ Fj ■ ■ - Fx, where the interval [0, t g ] is di- 



vided into N slices of duration At. One finds by eqn. [6] 



-At(iH(Ei(j))+r(Ei(j))) 



(8) 



with Ex(J) being the control amplitude in the 
j th time slice. The gradient of the fidelity 
can be computed in closed form as 



Retr<|Fj J F N ---F j+ i At 



a(ffl(E 1 q))+r(E 1 (j))). 



9Ei{j) 



■Fx 



We aim at optimizing the evolution of the qubit alone. 
Therefore, the TLF is traced out at the end of the full 
time evolution F(t g ) retaining all degrees of freedom in 
the intermediate steps Fj. One has schematically 



F li {t g )[]=tx^{F{t g )[.®p^]} 



(9) 



where the resulting map F R acts on the space of qubit 
density matrices alone. Here, we assume standard factor- 
ized initial conditions with the fluctuator in equilibrium, 
Ptlf- We use F R in fidelity eq. ^ for the optimization. 

Results and their discussion.— In the above model we 
focus on optimizing controls for a Z-gate, exp(— i%a z ). 
This is a paradigmatic case demonstrating the virtues of 
our technique: (1) an error rate up to approximately one 
order of magnitude lower than the current optimal work- 
ing point strategies; (2) the obtainable fidelities reach the 
Tx limit of the relaxation model; (3) the optimized con- 
trols relate to optimal times via self-refocussing effects — 
thus showing how open systems GRAPE-derived con- 
trols provide physical insight under structured relaxation 
models. Similar findings can be expected beyond one 
qubit gates. 

An overview of the accessible gate performance as a 
function of the duration t g of the gate is given in fig. [TJ 
(top). We notice that excellent gate performance can be 
achieved for pulse time t g > 7r/A. This corresponds to 
the static Aa x inducing at least a full loop around the 
Bloch sphere. Indeed, for the pulse at t g = 3. 375/ A 
the evolution consists of a quarter z-rotation, a full loop 
around x, and the second quarter of the z-rotation lead- 
ing to the total half rotation around z necessary for the 
Z-gate, see Bloch sphere in fig. [2] (right) for a particular 
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FIG. 1: Top: Gate error versus pulse time t g for optimal Z- 
gate pulses in the presence of a non-Markovian environment 
with dissipation strength k. A periodic sequence of minima 
at around t n = nir/A, where n > 1, is obtained. Middle: 
The gate error of optimized pulses approaches a limit set by 
Ti and 2Ti, as shown with k = 0.005. Bottom: Optimized 
pulses reduce the error rate by approximately one order of 
magnitude compared to Rabi pulses for n — 0.005. Pulses 
starting from zero bias and with realistic rise times (penalty) 
require only a small additional gate time. In all figures the 
system parameters are E2 = 0.1 A, A = 0.1 A and T = 0.2A. 



initial state. At shorter times the pulses cannot use the 
physical resource provided by the internal evolution to 
refocus the qubit. 

At longer times the attainable gate performance mildly 
deteriorates, depending on the value of k. This indicates 
that the optimal pulses are essentially limited by Ti pro- 
cesses at the optimal working point. We compare the 
performance to 1 — e~ tr ^ Tl with Ti obtained at E\ = 0. 
The optimized pulses are able to beat this limit which in- 
dicates that T\[E\ — 0) is a lower bound for the effective 
Ti, see fig. [T] (middle panel). For clarity, we also compare 
to l/ e -*/ T 2.""» with T 2>min = 2Ti. 

Optimizing the qubit with decoherence can lead to an 
improvement of around one order of magnitude over con- 
ventional Rabi pulses, see fig. [I] (bottom). The source 
of error of the Rabi pulses, can be seen as the fast 
oscillating counter-rotating component in the rotating 
frame, which is significant for these short, high-amplitude 
pulses, which only consist of a few carrier periods. Only 
at certain times t n ~ mr/A,n > 1, Rabi pulses perform 
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FIG. 2: Left: Rabi (- -) and optimized pulse ( — ) close to the 
optimal time at t g = 3. 375/ A and k = 0.05. Inset shows the 
evolution of initial state p = \+}{+\, ±) = (|0) ± \l))/V2 on 
the Bloch sphere under these pulses. Right: Comparison of 
pulse shape and qubit entropy between Rabi (- -), optimized 

( — ), and penalty-constrained pulse ( ) at t g = 5.0/A and 

k = 0. E 2 = 0.1A, A = 0.1A, T = 0.2A in all panels. 



well. At these t n we have n rotations around x due to the 
static field, refocussing the qubit. Additional optimiza- 
tion gives moderate results: in the case n = 2 the errors 
are (l-0) Rab i = 1.5-HT 3 and (1-0) G rape = 0.9-1CT 3 . 
At other times, the pulse optimization redirects the A 
drift to refocus the TLF-field and keeps the control at the 
optimal point E\ — as often as possible. As a result, 
the narrow and deep minima of the error for the regular 
Rabi pulse become shallow and broad using GRAPE. 

Steep rises and offsets of pulses as in fig. [2] can be 
smoothened by adding a penalty function to the fi- 
delity = 0- /J 9 a(t)Ef(t)dt. Although m princi- 
ple the penalty can be increased with the iterations 
in the algorithm, the simple penalty strength a(t) — 

2 — tanh (j^j + tanh was sufficient to avoid 
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offset and short rise times, see fig. [2] (right). Here, the 
overall penalty and the characteristic rise-time parameter 
are chosen to be a = 2A and t — 0.02/A, respectively. 
Fig. [T] shows that smooth pulses close to experimentally 
realistic settings come at a modest price of 0.5/A in addi- 
tional gate time. After t g » 3.75/A the same gate errors 
as without the penalty function are obtained. 

We now analyze the dependence on the bath parame- 
ters. In fig.[3]we can identify a nonmonotonic dependence 
of the error of the optimized pulse on 7. At low 7, we can 
approximate (1 — 0) =0 + 67 (left inset). First, the lin- 
ear growth with 7 accounts for the increasing probability 
for the TLF to flip at a random time during the evolu- 
tion. This is reflected in the power spectrum of the RTN 
for low 7, where essentially S(oj 3> 7) — A 2 j/lu 2 , eq. [5] 
Note that for time-independent E\ and very long evolu- 
tion, T 2 would be dominated by 5(0) in the Markovian 
limit. However, at very fast manipulations through an 
external field, the environment is only sampled at higher 



frequencies, 1 /t g 



^> 7. Second, it turns out to be 



that a ~ 0. This reflects that, for a static TLF (7 = 0) 
the GRAPE pulse fully compensates for the unknown ini- 
tial state of the TLF. At the end of the pulse, the entropy 
of the fluctuator is completely pushed out of the qubit, 
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lations, embodied in the Hamiltonian, turn out to be a 
crucial ingredient. 

We gratefully acknowledge support by NSERC discov- 
ery grants, by the EU in the projects EuroSQIP and QAP 
as well as by the DFG through SFB 631. 
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FIG. 3: Gate error versus TLF rate 7 for various temperatures 
for an optimized pulse with t g — 5.0/A. The left inset is a 
magnification of the low-7 P ar t °f the main plot and reveals 
the linear behaviour. The right inset shows the maximum of 
the curves of the main plot versus temperature. (E2 = 0.1 A 
and A = 0.1 A) 



see fig. [2] (right). 

On the other hand, for a high flipping rate 7, the 
physics of motional narrowing sets in. This limits the 
low frequency noise and hence the pure dephasing to 
S(u> <C 7) = A 2 /7 which vanishes for 7—5-00. In- 
deed, the high-7 part of the error can be fitted by a law 
c+ d/j. The finite limiting value c captures the residual 
decoherence which occurs even though the RTN model, 
eq. suggests absence of noise. Consequently, there 
is a 7 max at which the error is maximum. Remarkably, 
7 max ~ 0.32A ~ A/tt independent of any other param- 
eters such as temperature and pulse length: The attain- 
able performance is worst if the TLF flips once per free 
rotation around x. 

The maximal amplitude of the error (1 — 4>)("f m ax) as 
well as most other fit parameters show a non-monotonic 
temperature dependence, see fig. [3] (right inset), ft is ex- 
ponentially suppressed at low T and saturates to a finite 
high-temperature limit. The intermediate behaviour can 
be related to the equilibrium susceptibility of the TLF, 
which is maximum at T ~ 1Ei . The more responsive the 
TLF is to perturbations of its z-field, the stronger it will 
influence the qubit. 

Conclusion.— We have investigated an important 
model for decoherence in solid-state systems, a qubit cou- 
pled to a two-level fluctuator which itself is coupled to a 
heat bath. Our study is the first to exploit the explicit dy- 
namics of a complex non-Markovian environment in op- 
timal control of open systems for implementing quantum 
gates. The pulses include offset and rise time limitations 
of experimental settings. For a wide range of parameters, 
we have identified self-refocusing effects, which are usu- 
ally only visible at specific optimal pulse durations but 
can now be achieved more robustly. We have shown that 
both for fast and slow flipping of the TLF high-fidelity 
control can be achieved. The full qubit-fluctuator corre- 
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